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ABSTRACT 






Context. This study is part of a large project to study the physics of accretion and molecular outflows towards a selected 
sample of high-mass star-forming regions that show evidence of infall and rotation from previous studies. 
Aims. We wish to make a thorough study at high-angular resolution of the structure and kinematics of the HMCs and 
corresponding molecular outflows in the high-mass star-forming region G24. 78-1-0. 08. 
Q , Methods. We carried out SMA and IRAM PdBI observations at 1.3 and 1.4 mm, respectively, of dust and of typical 

$_( ■ high-density and molecular outflow tracers with resolutions of < 1". Complementary IRAM 30-m ^^CO and ^''CO 

"tx ' observations were carried out to recover the short spacing information of the molecular outflows. 

j^ I Results. The millimeter continuum emission towards cores G24 Al and A2 has been resolved into 3 and 2 cores, 

respectively, and named Al, Alb, Ale, A2, and A2b. All these cores are aligned in a southeast-northwest direction 
coincident with that of the molecular outflows detected in the region, which suggests a preferential direction for star 
formation in this region. The masses of the cores range from 7 to 22 Mq, and the rotational temperatures from 128 to 
^ I 180 K. The high-density tracers have revealed the existence of 2 velocity components towards Al, one of them peaks 

\l , close to the position of the millimeter continuum peak and of the HC Hll region, and is associated with the velocity 

gradient seen in CH3CN towards this core, while the other one peaks southwest of core Al and is not associated with 
C^ ' any millimeter continuum emission peak. The position-velocity plots along outflow A and the ^^CO (2-1) averaged 

^^ ' blueshifted and redshifted emission indicate that this outflow is driven by core A2. Core Al apparently is not driving 

f^ , any outflow. The knotty appearance of the highly coUimated outflow C and the ^'^CO position- velocity plot suggest an 

^— «, ■ episodic outflow, where the knots are made of swept-up ambient gas. 

Key words. ISM: individual G24. 784-0. 98 - ISM: molecules - ISM: jets and outflows - stars: formation 
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1. Introduction Keto & Zhang [^UTU]) . On the other hand, despite the fact 

5, , that theory predicts the existence of both large rotating 

^ . Molecular outflows, infall and rotation are ubiquitous phe- structures (e.g. Peters et al. MM and accretion disks (e.g. 

nomena duruig the earliest stages of the formation of stars Krumholz et al. mm around massive 0-type stars, such 

of all masses and luminosities. During the last few years, jisks have been elusive to date. One flnds, instead, huge 

much effort has been devoted to studying and describ- (_q^ p^)^ massive (a few 100 Mq) rotating toroids, which 

mg the physical properties of embedded massive protostars ^^^ ^^^ Keplerian and very likely non-equilibrium struc- 

and their molecular outflows (e.g. , Shepherd feChurchwell ^^^^^ (Beltran et al. mn\i . Finally, infall onto young 0-type 

1996a,b; ZhMLg et al. i2001j Beuther et al. [2DD2 Lopez- g^^^g ^^^ ^een measured through high-angular resolution 

Sepulcre[2DDi). The presence of rotating structures around observations of lines showing redshifted self-absorption or 

massive young stars has also been asserted thanks to the absorption against a background bright continuum source 

detection of velocity gradients perpendicular to the axis (j^eto et allMIl Zhang et al.[l998, Behran et aL [MM|l . or 

of the molecular outflows m the cores (see Cesarom [20021 by measuring the velocity field of the ionized gas surround- 

and reference therein). These studies have revealed two j^^g ^ newly formed massive star (see Keto [Mi- 
different kinds of structures surrounding B- and 0-type 

stars. On the one hand, B-type stars are surrounded by G24.78+0.08 (hereafter G24) is one of the most studied 
centrifugally supported disks, as seen in the prototypical massive star-forming regions. The simultaneous presence of 
B-type protostar IRAS 20126-1-4104 (Cesaroni et al. 120051 the three elements fundamental for star formation (outflow, 
infall, and rotation) has been reported in one of the mas- 
Send offprint requests to: M. T. Beltran, e-maU: sive objects embedded in this clump. This region, which 
mbeltrELii@arcetri.astro.it is located at a distance of 7.7 kpc, has been studied in 
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Table 1. Parameters of the maps 
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" I'he synthesized CLEAN ed beams for maps with the KOBUSl" parameter ot Briggs (I1995II set equal to 0. 




' For the molecular line observations the la noise is 


per channel. 










" The synthesized beam for the SMA VEX configuration map is 0'.'55 x 0'.'40 at P.A.=52° 


, and the rms noise is 


3.2 mJy beam~^. 



The frequency is that of the K=Q component. 

The synthesized CLEANed beams for maps using natural weighting. 



great detail by our group (Codella et al. 119971 Furuya et 
al. lM?^ Cesaroni et al .lMJgl Beltran et al .BUM Beltran et 
al. [20051 B eltran et al.lMISl Beltran et al. lMTTI Moscadelli 
et al.HnnZl Vig et al. [M?5)) . G24 contains a cluster of Young 
Stellar Objects (YSOs), three of which have rotating toroids 
associated with them (Beltran et al. l2004p . These rotating 
toroids seem to be associated with two compact molecu- 
lar outflows, first mapped in CO by Furuya et al. (j2002p . 
In one of these rotating cores, G24 Al, the presence of an 
early-type star is witnessed by an embedded hypercompact 
(HC) Hil region, first detected at 1.3 cm by Codella et 
al. (|1997p and later studied at very high-angular resolution 
at 1.3 and 0.7 cm by Beltran et al. (|2007p . The properties 
of the ionized gas have also been studied through contin- 
uum emission and recombination lines by Galvan-Madrid 
et al. (pnUS)) and Longmore et al. (^005)) . This HC Hii re- 
gion, which has a diameter of ^0'.'2 (^1,500 AU), is located 
at the geometrical center of the toroid. The free-free con- 
tinuum spectrum resembles that of a classical (Stromgren) 
Hii region around a zero-age main sequence star of spectral 
type 09.5, corresponding to a mass of ~ 20 Mq. Beltran 
et al. (|2006p mapped the core G24 Al in NH3 and detected 
redshifted absorption, indicating that the toroid is under- 
going infall towards the HC Hii region. 



We have carried out new 1.3 mm interferometric obser- 
vations of this region at sub-arcsecond resolution with the 



Submilhmeter Array (SMAfl (Ho et al. [MH)) of several 
outflow and high-density tracers. The goals of this study 
were to better characterize the structure and kinematics 
of the Hot Molecular Cores (HMCs) and molecular out- 
flows in this region, and in particular, to identify the pow- 
ering source of the outflow observed towards G24 Al. In 
order to be sensitive to extended structures filtered out by 
the interferometer, we have also mapped the region over 
2' with the IRAM 30-m telescope). In addition, we used 
the enhanced capabilities of the IRAM Plateau de Bure 
(PdBI) to observe this region in ^^CHsCN (12-11) at very 
high-angular resolution to investigate whether the massive 
toroids were hiding true (Keplerian) circumstellar disks in 
their interiors. Unfortunately, the angular resolution of the 
observations has not been enough to detect circumstellar 
disks embedded in the massive toroids, but it has allowed 
us to study the velocity structure of the cores with un- 
precedented resolution. The observational details are given 
in Sect. 2, while the results are illustrated in Sect. 3 and 
analyzed and discussed in Sect. 4. Finally, the conclusions 
are drawn in Sect. 5. 



* The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded 
by the Smithsonian Institution and the Academia Sinica. 

2 IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 
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Fig. 1. (Left panel) SMA map of the 1.3 mm continuum, emission obtained combining both compact and very extended 
configurations. Contour levels are 3, 6, 9, 12, 15, 35, and 60 times a, where la is 3 mJy beam~^. The crosses mark the 
positions of the sources identified by Furuya et al. (|2002) and Beltran et al. (|2004p . The triangle marks the position of 
the HC Hii region in G24 Al (Beltran et al. I2007p . The SMA synthesized beam is shown in the lower left corner. {Right 
panel) Very extended configuration only SMA map at 1.3 mm (contours) overlaid on the PdBI map at 1.4 mm (grayscale) 
towards G24 Al, A2 and D. Contour and grayscale levels are 3, 6, 9, 12, 15, and 20 times a, where la is 3.2 mJy beam~^. 
The SMA synthesized beam is shown in the lower left corner and the PdBI one in the lower right corner. 



2. Observations 

2.1. Submillimeter Array 

G24 was observed with the SMA using two different array 
configurations. Compact-array observations were taken on 
July 17, 2007, and covered baselines with lengths between 
7 and 100 kA (sampling spatial structures in the range of 
29'.'5 to 2'.'1). Very Extended (VEX) configuration data were 
made on August 2, 2008, with baseline lengths from 23 to 
391 kA (9'.'0 to 0'.'5). The SMA had two spectral sidebands 
- both 2 GHz wide - separated by 10 GHz. For both obser- 
vations, the two sidebands covered the frequency ranges of 
219.3-221.3 and 229.3-231.3 GHz with a uniform spectral 
resolution of '^0.5 km s^^. 

The phase reference center of the observations was set 
to the position a(J2000)=18^ 36"^ 12!565, S{3200Q)=-Q7" 
12' 10'.'90. Absolute fiux calibration was derived from ob- 
servations of Callisto. The bandpass of the receiver was cal- 
ibrated by observations of the quasar 3C454.3. Amplitude 
and phase calibrations were achieved by monitoring 1733- 
130. The flux densities estimated for 1733—130 are 1.71 and 
2.29 Jy for the compact and the VEX configuration, respec- 
tively. We estimate the fiux-scale uncertainty to be better 
than 15%. The visibilities were calibrated with the IDL su- 
perset MlBp. Further imaging and analysis was done with 
MIRIAD (Sault et al. 1995) and GILDA^. The contin- 
uum was constructed in the ('Wj^^j-domain from the line- free 
channels. Continuum and channel maps were created com- 
bining the data of both compact and VEX configurations 
with the ROBUST parameter of Briggs (|1995p set equal to 



^ The MIR cookbook by Charlie Qi can be found at http:// 
cfa-www.harvard.edu/'^cqi/mircook.html 



The 



GILDAS 



package 



available 



at 



http://www.iram.fr/IRAMFR/GILDAS 



0, except for the ^^CO, i^co, Ci^O, and SO channel maps 
that were created using natural weighting. Details of the 
synthesized CLEANed beam, spectral resolution, and rms 
noise of the maps obtained in a number of lines are given 
in Table m 

2.2. IRAM Plateau de Bure Interferometer 

PdBI observations at 214.3 GHz in the B configuration were 
carried out on March 9 and 11, 2007. The correlator was 
set up to observe the ^^CHsCN (12-11) line emission at 
216 GHz with both polarizations. The units of the correla- 
tor were placed in such a way that a frequency range free 
of lines could be used to measure the continuum flux. A 
unit of 80 MHz of bandwidth was set to cover the emis- 
sion of "CH3CN (12-11) K =0 to 4, with a spectral res- 
olution of 0.156 MHz or 0.22 km s~^. A second unit of 
320 MHz of bandwidth was set to cover the emission of 
"CH3CN (12-11) K =0 to 7, with a spectral resolution 
of 2.5 MHz or ~3.5 km s^^. The two remaining units of 
320 MHz of bandwidth were placed to cover the rest of the 
band. 

The phase center used was a(J2000)=18h 36'" 12!661, 
(5(J2000)=-07" 12' 10'.'15. The bandpass of the receivers 
was calibrated by observations of the quasar 3C273. 
Amplitude and phase calibrations were achieved by moni- 
toring 1741—038, whose flux densities were determined rela- 
tive to MWC349 or 1749-1-096. The flux densities estimated 
for 1741-038 are 1.18 and 1.12 Jy for March 9 and 11, 
respectively. The uncertainty in the amplitude calibration 
is estimated to be ~20%. The data were calibrated and 
analyzed with the GILDAS software package. The contin- 
uum maps were created from the line free channels. We 
subtracted the continuum from the line emission directly 
in the (^M^wj-domain. Continuum and channel maps were 
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Table 2. Parameters of the cores 





Position 


peak*^ 
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Source Diameter"''' 


Trot 
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° The positions and diameters have been estimated from the very extended configuration SMA map. 

' Deconvolved average diameter of the 50% contour level. 

'^ From the rotational diagram method (see Sect. 14. 1|) . For those cores for which an estimate of Trot is not possible with this 

method, a range of temperatures is given. 

'^ Mass estimated using a dust temperature equal to Trot and a dust opacity of 0.8 cm^ g~^ at 1.3 mm (see Sect. 14. l]) . 

"^ The dust continuum emission associated with core Al and used to estimate the mass is 102 mjy (see Sect. 13. l|) . 

' The dust continuum emission associated with cores Al+Alb+Alc and used to estimate the mass is 780 mJy (see Sect. l3.l"|l . 



created using natural weighting. Details of the synthesized 
CLEANed beam, spectral resolution, and rms noise of maps 
are given in Table [T] 

2.3. IRAM 30-m telescope 

Supplementary short-spacing information to complement 
the interferometric data was obtained with the IRAM 30-m 
telescope. The observations were carried out on November 
11th, 2006 with the HERA receiver using the On-The- 
Fly mode. The receiver was tuned to 230.538 GHz for 
i2CO(2-l) (HERAl) and 220.399 GHz for ^^CO (2-1) 
(HERA2). Both lines were observed simultaneously and 
covered with the VESPA autocorrelator with 0.1 km s~^ 
spectral resolution, corresponding to a channel spacing of 
~0.07 MHz. The maps sizes were 2' x 2' and were centered 
at a(J2000)=18'^36'"12!6 and (5(J2000)=-07°12'10'.'9. The 
sampling interval was 2" and the region was alternatively 
scanned in the north-south direction and the east-west di- 
rection in order to reduce effects caused by the scanning 
process. 

The spectra were reduced using CLASS of the GILDAS 
software package and channel maps with 1 km s~^ resolu- 
tion were created. The ^^CO data have a beam size of 10'.'4 
and the "CO of 10'.'9. The IRAM 30-m data sets were com- 
bined with the SMA data sets using the UVMODEL task 
of the MIRIAD package. The synthesized beam of the com- 
bined SMA-hlRAM 30-m data is 3'.'37 x 1'.'96 at P.A.=70° 
for i^CO and 3'.'57 x 2'.'02 at P.A.=71° for ^^CO 



3. Results 

3.1. Continuum emission 

Figure [T] shows the SMA map of the 1.3 mm continuum 
emission towards G24. The inset shows the map obtained 
only from the very extended configuration at 1.3 mm to- 
wards G24 Al, A2, and D overlaid on the PdBI continuum 



map at 1.4 mm. The positions, fluxes, and deconvolved sizes 
of the cores measured as the average diameter of the 50% 
contour from the SMA map are given in Table [2l 

As seen in Fig. [H all but core B have been detected at 
1.3 mm. The emission towards cores Al and A2 has been 
resolved into additional cores with the higher angular reso- 
lution provided by the very extended configuration observa- 
tions (Fig.[T]). Core Al has been resolved into 3 cores named 
Al, Alb, and Ale, whereas core A2 has been resolved into 
two named A2 and A2b. The differences between the posi- 
tion and geometry of cores Ale and A2b in the SMA VEX 
configuration map and the PdBI one are probably due to 
the slightly extended structure of the sources and the differ- 
ent (u,v)-coverage of the observations. The cores are aligned 
in a southeast-northwest direction coincident with that of 
the molecular outflow associated with cores Al and/or A2 
(see Figs. [5] and O. The core associated with the HC Hii re- 
gion is Al. To derive a rough estimate of the free-free contri- 
bution of the embedded HC Hll region at millimeter wave- 
lengths, we extrapolated the optically thin free-free 7 mm 
continuum emission {Si, oc i/~'^'^) measured by Beltran et 
al. (|2007p . The integrated flux density at 7 mm is 101 mJy, 
and the expected free-free emission at 1.3 mm is ~85 mJy. 
Therefore, the dust continuum emission associated with Al 
is 102 mJy (see tabled]). Note that the continuum peak in- 
tensity of core Al is ~86 mJy/beam, which is very similar 
to the expected free-free emission. This means that there is 
almost no dust emission within the HC Hii region. 

We have also detected, although at a 3a noise level, core 
D at 1.3 mm with the SMA and 1.4 mm with the PdBI (see 
Fig. [T|). This core was first detected at 2.6 mm with the 
PdBI and at 2 mm with the Nobeyama Millimeter Array 
(NMA) by Furuya et al. (|2002)) . Core D was not detected at 
3.3 and 1.4 mm with the PdBI by Beltran et al. (2001]) be- 
cause the expected flux at those wavelengths fell below the 
3a noise level of the observations. In our case, the detection 
of core D is very weak, the measured SMA flux density is 
14 mJy with the VEX configuration, and ^^^16 mJy with 
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Fig. 3. ^^CHgCN (12-11) spectrum obtained with the PdBI 
by integrating the emission inside the Scr contour level area 
towards G24 Al and A2. Different K numbers are marked 
with dashed lines in the upper part of the spectrum. The 
spectral resolution has been degraded to 0.5 km s~^. 



spectral resolution of 0.25 km s'^ The "CH3CN (12-11) 
K — h and 6 components, which have been observed with 
a spectral resolution of 3.5 km s~^, have also been detected 
towards both cores G24 Al and A2. Maps of the CH3CN 
(12-11) emission averaged under the K = 2 and 6 compo- 
nents, and of the CHg^CN (12-11) and i^CHgCN (12-11) 
emission averaged under the K = 2 component towards 
G24 are shown in Fig. |4l 

As seen in Fig. 21 the CH3CN and isotopologues emis- 
sion is clearly detected and resolved towards both cores Al 
and A2. The emission peaks at the same position as the 
millimeter continuum. 



3.3. Other high-density tracers 

In addition to CH3CN and isotopologues, most of the 
species detected towards G24 are high-density tracers. 
Figure [5] shows the emission averaged in the velocity in- 
terval (105-115) km s~^ for the high-density tracers. The 
molecular outflow tracers are discussed in the next section. 
As seen in this figure, all the species have been clearly de- 
tected towards cores Al and A2. All the species peak to- 
wards the millimeter continuum peaks. On the other hand, 
only CH3OH, 13CS, HNCO, OCS, and Ha^^CO have been 
detected towards core C, although in most cases the emis- 
sion is so weak that it is not possible to distinguish any 
clear emission peak. 



both configurations (compact plus VEX). This latter value 
is consistent with the extrapolation of the fit to the spec- 
tral energy distribution by Cesaroni et al. (|2003p . assum- 
ing that the source size is similar to the synthesized beam 
size of the SMA observations as indicated by the contin- 
uum maps. The fact that this source is almost completely 
resolved out with a resolution of '^0'.'84 suggests that the 
source has two components: a compact core, which is what 
we have mapped at 1.3 and 1.4 mm, plus an extended halo 
whose emission is filtered out by our observations. 

3.2. CH3CN and isotopologues 

Figure [5] shows the spectra of the upper and lower SMA 
sideband around 230 and 220 GHz, respectively, obtained 
by integrating the emission inside the Scr contour level area 
towards G24 A1+A2, and G24 C. Figure [3] shows the PdBI 
spectrum around 214.3 GHz obtained by integrating the 
emission inside the 3a contour level area towards G24 Al 
and A2. In Fig. [51 only the lines that have been analyzed are 
labeled. These are the lines with a good signal-to-noise and 
that could be unambiguously identified, i.e., that are not 
blended with other lines. As seen in these figures, altogether 
we detected 11 species (CO, HNCO, Ha^CO, SO, CH3OH, 
HCOOCH3, CH2CO, CH3CN, C2H5CN, OCS, "CS), in- 
cluding 2 CO isotopologues (Ci^O and ^^CO) and 2 CH3CN 
isotopologues (CH3I3CN and 13CH3CN). 

Regarding the methyl cyanide emission, several K- 
components of CH3CN (12-11) (up to K=8 for Al and A2) 
are clearly detected towards G24 Al, A2, and C (Fig. EJ. 
Several if -components of CH3I3CN and 13CH3CN (12-11) 
(up to K=6) are also clearly detected towards G24 Al and 
A2 (Figs. H and EJ. Figure H shows only the if = to 
4 components, which are those observed with the highest 



3.4. Molecular outflow tracers 

The ^^CO (1-0) emission towards G24 has been previously 
studied by Furuya et al. (|2002|) . These authors have discov- 
ered two bipolar molecular outfiows in the region, oriented 
in the same direction, northwest-southeast: one associated 
with cores Al and/or A2 and named outflow A, and the 
other one associated with core C and named outflow C. 
The cloud systemic velocity, Vlsr, is 111 km s~^ (Furuya 
et al. I2002p . As seen in Fig. [2l we have not only clearly 
detected ^^CO (2-1) towards both cores, but also the CO 
isotopologues ^^CO (2-1) and Ci^O (2-1), and SO (65- 
54). Figures [6] and [7] show the channel maps of the ^^CO 
and "CO (2-1) SMA+IRAM 30-m combined emission. 
Note that the extended emission seen in the -'^^CO (2-1) 
SMA-I-IRAM 30-m combined channel maps for velocities 
between 82 and 86 km s~^ is due to a gas component along 
the line-of-sight. This component is not seen in ^^CO. As 
seen in the channel maps, the outflow A is clearly seen in 
both tracers, while the outflow C is clearly visible in ^^CO 
and only marginally in ^■^CO. 

Figure [5] shows the ^^CO, "CO, C^^O (2-1), and SO 
high-velocity blueshifted and redshifted averaged emission 
observed with the SMA. The maps shown in this flgure have 
been reconstructed with a circular beam of 2" to enhance 
the extended emission. As seen in this figure, the outfiow 
A, which has a P.A.=140°, is clearly seen in all lines. On 
the other hand, the outflow C, which has a P.A.=145°, is 
clearly visible in ^^CO and ^^CO and just barely visible 
in C^^O (2-1). The high- velocity SO (65-64) emission is 
hardly seen at redshifted velocities towards core C. 

As seen in Fig. \E\ the molecular outflow A shows a 
very prominent and poorly collimated blueshifted lobe and 
a more collimated redshifted lobe. This is especially visi- 
ble in the higher angular resolution ^^CO map visible in 
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Fig. ini The redshifted lobe shows an extension towards the 
northwest and the blueshifted one towards the southeast, 
especially in ^-^CO. The NW redshifted emission, detected 
at velocities close to the Vlsr, is still clearly visible for ve- 
locities higher than 40 km s~^ with respect to the Vlsr- 
The redshifted NW emission can be seen in the position- 
velocity plots (Fig. ^U\j from about —6" to —4", and the 
blueshifted SE emission from about 1" to 3". A possible ex- 
planation for this blueshifted (redshifted) emission towards 
the redshifted (blueshifted) lobe of outflow A could be that 
outflow A is close to the plane of the sky. Alternatively, it 
could indicate the presence of a second outflow in this area. 

The high-angular resolution maps (Fig.|9]) and position- 
velocity plots (Fig. [TU| seem to suggest that the core pow- 
ering this molecular outflow is A2, because it is located 
closer to the geometrical center of the outflow. However, 
with these observations we cannot exclude the possibility 
of the existence of more than one outflow (see Sect. 14.3.11 
for a more detailed discussion). 

The molecular outflow C is highly collimated. This can 
be seen in the ^^CO maps, either the 2" resolution map or 
the high-angular resolution one. As clearly seen in Figs. 13 
[21 andlini the blueshifted lobe is much smaller than the red- 
shifted one. This is particularly evident in the blueshifted 
emission maps averaged for high outflow velocities (not 
shown here), which we considered as those higher or lower 
than 36 km s~^ with respect to the Vlsr- For this high 
outflow-velocity range, the size of the blueshifted lobe is 
~3 times smaller than that of the redshifted lobe. Taking 
into account the fact that the molecular outflow consists of 
ambient gas that is swept up from a position close to where 
it is observed, a possible explanation for the morphology 
of the outflow C could be that the blueshifted lobe is close 
to the edge of the molecular cloud. That is, the blueshifted 
wind passes into a region where there is no, or little, molec- 
ular material to be swept up. 

Another feature of the outflow C is the redshifted emis- 
sion seen southeastern of core C. This redshifted emission 
is like a finger pointing towards the east, very prominent 
for outfiow velocities close to the Vlsr, but still visible for 
outflow velocities higher than 40 km s""'^ with respect to 
the Vlsr- This feature has a length of -8'.'4 (-0.31 pc). 
A possible explanation for this SE redshifted feature could 
be that the molecular outflow southwards of core C is in- 
teracting with the ambient cloud close to the position of 
the core. As a result of this interaction, the molecular out- 
flow would be deflected towards the east and the emission 
would shift from blueshifted to redshifted. A similar sce- 
nario has been explained by Beltran et al. (|2002p for the 
intermediate-mass outflow powered by IRAS 21391-1-5802 
in terms of the shocked cloudlet model scenario, where an 
inverted bow shock is produced in the windward side of a 
dense ambient clump fSchwartz ll978p . This scenario could 
also explain the smaller size of the blueshifted lobe with 
respect to the redshifted one, because part of the material 
that would be entrained by the molecular outflow south- 
wards of core C, would be deflected and redshifted. We 
have searched for a possible high-density southern clump 
responsible of the deflection of the outflow and have found 
a clump in CH3OH (8i,8-7o,7)E at a position a(J2000)=18'^ 
36™ 13'.'09, (5(J2000)=-07'^'l2™ T.7. However, the clump is 
only seen for redshifted velocities, ~2-3 km s"""^ with re- 
spect to the Vlsr, while one would expect the clump to be 
blueshifted like the southern lobe before deflection. Another 



possible explanation for the SE redshifted emission could be 
the presence of a second outflow, powered by a secondary 
source located close to core C but not resolved by our ob- 
servations, but the blueshifted lobe would be missing. 

Taking the blueshifted and redshifted emission into ac- 
count, the total extent of the outflows is — 13'.'5 (—0.5 pc) 
for outflow A and '-15'.'2 (-0.57 pc) for outflow C. 



4. Analysis and Discussion 

4.1. Temperature and mass estimates 

CH3CN is a symmetric top molecule, so it can be used 
to estimate the temperature in the cores. The rotational 
temperature, Trot, and the total methyl cyanide column 
density, A^cHsCN, can be estimated by means of the rota- 
tion diagram method, which assumes that the molecular 
levels are populated according to LTE conditions at a sin- 
gle temperature Trot- In the high-density limit where level 
populations are thermalized, one expects that Trot = Tkin, 
the kinetic temperature. The CH3CN ground-level tran- 
sitions appear to be optically thick for all cores, as sug- 
gested by the ratio between the main species and isotopo- 
logues. Therefore, the Boltzmann plot was only performed 
using the ^"^CHsCN and vibrational excited CH3CN transi- 
tions (see Fig. ITT]) . The core Ale has not been detected in 
13CH3CN and only marginally in CH3I3CN. Therefore, the 
rotation diagram in Fig. [TT] has been calculated with only 
three measurements, one corresponding to CH3^'^CN K=2, 
and two corresponding to CH3CN vg, = 1, {K,l)—{3,~1) 
and {K, 1)={6,+1). The Trot estimated for this core is 146 K, 
and A^CHaCN is about 1x10^^ cm~^. The relative abun- 
dance [CH3CN]/[i3CH3CN] was estimated to be 46 follow- 
ing Wilson & Rood (il994|) for a galactocentric distance 
of 3.7 kpc. The temperatures and column densities ob- 
tained are consistent with those estimated by Beltran et 
al. (2005) from CH3"CN (12-11) and (5-4), and CH3CN 
(5-4) W8 = 1- 

The masses of the individual cores, which are given 
in Table [21 have been estimated from the 1.3 mm dust 
continuum emission observed with the SMA at the high- 
est angular resolution (VEX) assuming a dust opacity of 
— 0.8 cm^g~^ at 1.3 mm (Ossenkopf & Henning |1994p . a 
gas-to-dust ratio of 100, and a dust temperature equal to 
Trot- For the three cores for which an estimate of Trot is 
not possible with the rotation diagram, a range of masses 
is given assuming that the dust temperature of the cores 
ranges from 30 to 60 K. The lower value is the kinetic tem- 
perature of —30 K estimated from ammonia observations 
carried out with the Very Large Array (VLA) by Codella et 
al. (I1997P towards G24. However, taking into account that 
CH3CN has been detected towards core C (Fig. |31), we ex- 
pect that the temperature, at least for this core, is higher. In 
Tablej^lwe also give the masses estimated from the compact 
plus VEX configuration SMA map. In this case, the emis- 
sion of some individual cores cannot be resolved, and there- 
fore, the masses correspond to that of cores Al-|-Alb-|-Alc 
and A2+A2b. 



4.2. Dense gas kinematics 

Following Behran et al. ((20041 [20051 [20TTll . the velocity 
field of the cores has been studied thanks to the CH3CN 
and isotopologues emission. 
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4.2.1. Velocity gradients 

Figure [12] shows the map of the CH3CN Hne velocity for 
the G24 Al and A2 cores overlaid on the VEX contin- 
uum emission. The line velocity map has been obtained 
by simultaneously fitting the CH3CN K = 0,1,2,3, and 
4 components, assuming identical line widths and fixing 
their separations to the laboratory values, at each position 
where CH3CN is detected. The validity of the assumption 
of identical line widths has been checked by fitting sepa- 
rately different JsT-components. Beltran et al. (2005) cal- 
culated very small variations in the CH3CN (12-11) line 
widths from i^ = to 6 (see their Table 7). As seen in 
Fig. [121 we have clearly detected the same velocity gradients 
in the southwest-northeast direction observed by Beltran et 
al. (2004) towards G24 Al and A2. 

The continuum emission of core Al peaks at the cen- 
ter of the southern velocity gradient, which is very close to 
the position of the HC Hii region. The P.A. of this velocity 
gradient is about 50°. The velocity shift measured over an 
extent of '^13000 AU is ~3.4 km s~^, that is, a velocity gra- 
dient of about 54 km s~^ pc~^. The peak of the continuum 
emission of core A2 is close to the center of the northern 
velocity gradient, which has a P.A. of ~40°. The velocity 
shift measured over an extent of ^11500 AU is ^2.4 km s^^, 
that is, a velocity gradient of about 43 km s^^ pc^^. 

Unfortunately our data do not have enough sensitivity 
to allow us to properly map the velocity field towards the 
core C and better constrain the direction of the velocity 
gradient detected in CS (3-2) by Behran et al. p004| . 

4.2.2. Multiple components towards G24 Al 

The analysis of "CH3CN (12-11) has clearly revealed the 
existence of 2 velocity components towards core Al. This 
is clearly seen in Fig. [2| for K='i and 4. The second ve- 
locity component is less evident for K—1 because it is 
blended with K—0, and for K=2 because it is blended with 
some species typical of hot cores, such as dimethyl ether 
(CH3OCH3) and/or methyl formate (CH3OCHO). On the 
other hand, the K=Q transition is very close to the edge of 
the observed band, which makes it very difficult to distin- 
guish the second velocity component. The 2 velocity compo- 
nents can also be seen in CH3^'^CN and C2H5CN as shown 
in Fig. [121 where we show the spectra averaged over the 
core towards Al and A2. Note that the reason for showing 
CHa^CN (12-11) K=2 instead of K=Z is due to the fact 
that the latter one is blended with CH3CN (12-11) K=Q. 
As seen in Fig.[T31 the spectra towards Al clearly show 2 
components. Component 1 that peaks at about 106 km s^^ 
(107 km s^^ for C2H5CN), and component 2 that peaks 
at about 111 km s~^ (the systemic velocity). The spectra 
towards A2 do not show an apparent component peaking 
at a velocity of '^106 km s~^, although they do show a 
significant blueshifted wing. The bottom panels in Fig. [T51 
show the "'^^CH3CN (12-11) K — 3 emission averaged for 
a velocity interval (105, 107.5) kms~^, corresponding to 
component 1, and (108.75, 114.25) km s^^, corresponding 
to component 2. As seen in these maps, the 2 components 
peak at a different position. Component 1 peaks towards 
the southwest of the core Al continuum emission peak, 
and its emission is clearly elongated in the SW direction. 
On the other hand, the emission peak of component 2 coin- 
cides with that of the core Al continuum emission and with 



the position of the HC Hii region. For core A2, the peak 
of the emission averaged for the 2 velocity intervals is sim- 
ilar although not coincident, and in both cases is displaced 
from the position of the continuum emission peak. A SW 
elongation is also visible for component 1, but not for com- 
ponent 2. Hence, we cannot discard a priori the existence 
of 2 velocity components also for core A2. 

We have searched for the two different velocity com- 
ponents in the other high-density tracers observed with 
the SMA, and in particular in CH3CN and CH3CN U8=l. 
However, although the spectral resolution of the observa- 
tions is enough to disentangle the emission of the two com- 
ponents, what we have found is a significant blueshifted 
wing in some cases, but not two clearly separated com- 
ponents. We have reconstructed the PdBI ^^CHsCN K=3 
and the SMA CH3CN K=3 maps with the same synthe- 
sized beam of 0'.'89 x 0'.'0.59 P.A. =10°, i.e. a beam encom- 
passing the beams of the SMA and PdBI maps, and have 
recalculated the spectrum towards Al. The component 1 
was still clearly visible in -'^^CH3CN, even after degrading 
the spectral resolution to that of the SMA observations 
(0.6 km s"i), and in CHs^^cj^ Qn the other hand, the 
component 1 was still not evident in the CH3CN spectrum. 

We have studied the velocity field of the 2 components 
in ^'^CH3CN towards Al to see whether the presence of 
another velocity component at ~106 km s~^ could be re- 
sponsible for the SW-NE velocity gradient seen in CH3CN 
towards core Al (see Fig. [12]). In fact, a velocity gradient 
might also be mimicked by two distinct cores with different 
^LSR that are too close to be resolved by our observations. 
Figure [14] shows the first moment map of C2H5CN com- 
puted over a velocity interval from 108.5 to 114.5 km s~^, 
that is, over the range of the second velocity component. 
We have analyzed the velocity field of C2H5CN, because 
this is the high-density tracer with 2 components clearly 
distinguishable that shows the strongest emission towards 
G24 Al (see Fig. [13]) . The analysis indicates that compo- 
nent 2 shows a velocity gradient similar to the one observed 
in CH3CN, centered close to the position of the contin- 
uum emission peak and the HC Hll region, although with 
a velocity range slightly smaller (~109.9-111.4 km s~^). 
Therefore, we conclude that the velocity gradient seen to- 
wards core Al, first detected in CH3CN by Beltran et 
al. (|2004p . could only be marginally affected by the presence 
of a second core at a different Vlsr, and is produced by a 
structure rotating around the star powering the HC Hii re- 
gion. Hence, the velocity component at ^106 km s^^ would 
not be related to this gradient. The presence of a SW core 
is also discarded by the fact that a secondary continuum 
emission peak is not detected in that direction (Fig. [Q]) . 
One would expect that any embedded core should be bet- 
ter detected in continuum than in line emission. 



4.3. The molecular outflows 

Figures \8\ and [9] show the molecular outflows A and C 
mapped in ^^CO (2-1), ^^co (2-1), C^^O (2-1), and 
SO (65-54) towards G24. In Fig. [TU]we show the position- 
velocity plots along the outflows A and C. For the sake 
of completeness we show both the SMA and the combined 
SMA-flRAM 30-m data. The outflow C is hardly visible 
in i^CO (2-1) with the SMA at a resolution of ~0'.'7, and 
hence, for this outflow we show only the ^^CO (2-1) data. 
Note that the cuts were made along P.A. =—40° (outflow 
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Table 3. Properties of the molecular outflows"'' 



Outflow 


R"^ cos i 
(pc) 


Mout 

(Mo) 


Mout tan i 


P"^ sin i 

{Mq km s"i) 


(10"« erg) 


F"^ sin^ i/ cos i 
(M0kms"Vr"') 


tout cot i 

(yr) 


A 

C red finger-'^ 


0.50 
0.57 
0.31 


27.3 
7.3 
4.2 


1.5 X 10"'' 
3.0 X 10"" 
2.2 X 10"" 


245 
62 
33 


2.5 
0.6 
0.3 


1.4 X lO""" 

2.6 X 10"^ 

1.7 X 10"^ 


1.8 X 10" 
2.4 X 10" 

1.9 X 10" 



(a) The tabkle indicates the estimated values. In case of knowing the inchnation angle with respect to the plane of the sky, i, of 
the outflows, the parameters should be corrected accordingly. 

(b) The blueshifted emission has been integrated for the velocity range (94, 103) km s"^, and the redshifted one for (116, 
127) kms"\ 

(c) Total size of the lobes. 

(d) Momenta and kinetic energies are calculated relative to the cloud velocity. 

(e) The redshifted emission seen SE of core C has not been taken into account in the calculations. 

(f) Redshifted emission seen SE of core C. 



A) and P.A.=— 45° (outflow C) and all the plots have been 
obtained after averaging the emission along the direction 
perpendicular to the cut with the purpose of increasing the 
S/N of the plots. 

4.3.1. Which source is powering outflow A? 

One of the questions that still remains open in the study 
of G24 is which core, Al or A2 (or both), is driving out- 
flow A. The high-angular resolution maps (Fig. [S]) seem to 
indicate that the core powering this outflow is A2, because 
it is located closer to the geometrical center of the outflow. 
This is especially visible in the ^'^CO map, where outflow A 
shows a clear bipolarity centered in core A2. However, the 
most clear evidence that A2 is the powering source comes 
from the position- velocity plots along the outflow (Fig. [TU)) . 
In these plots the vertical solid and dashed lines indicate 
the position of core Al and A2, respectively. As can be 
seen in these plots, core A2 is located at the position where 
the emission, or the emission at the highest outflow ve- 
locity, changes from being blueshifted to redshifted. This is 
more evident in the ^"^CO SMA plots, because no redshifted 
emission is detected at any position before that of core A2. 
Regarding core Al, the position- velocity plots seem to in- 
dicate that this core is not powering any outflow. However, 
as shown in Fig. 1^1 the CO emission is very complex to- 
wards cores Al and A2, with redshifted emission towards 
the blueshifted lobe and the opposite. Hence, we cannot to- 
tally discard the possibility that core Al could be powering 
an additional outflow in the region. In any case, this pos- 
sible outflow would not be the most prominent one, which, 
as already mentioned, is driven by core A2. 

4.3.2. Physical parameters 

Table [3] gives the parameters of outflows A and C, as well 
as those of the redshifted emission feature seen SE of core 
C. The emission of this feature has not been taken into ac- 
count when calculating the parameters of outflow C. The 
parameters have not been corrected for inclination angle, i, 
of the flow with respect to the plane of the sky. However, we 
have indicated in Table [3] the correction for the inclination 
to be applied for each outflow parameter. The size of the 
lobes, mass Mout, outflow mass loss rate Mout, momentum 
P, kinetic energy E, momentum rate in the outflows F, and 
dynamical timescale tout were derived from the ^'^CO emis- 
sion for the velocity ranges indicated in Table |21 The ^'^CO 



channel maps used were those reconstructed with a 2" cir- 
cular beam. The dynamical timescale of the blueshifted and 
redshifted lobes was estimated as iout = R/Vout, where Vout 
is the difference in absolute value between the maximum 
blueshifted or redshifted velocity and the systemic velocity, 
Vlsr=111 km s~^. The iout of the outflow is the maximum 
dynamical timescale of the two lobes. The parameters R, 
Mout, P, and E were calculated for the blueshifted and red- 
shifted lobes separately, and then added to obtain the total 
value. The ["C0]/[H2] abundance rati o of 2 2 x 10"*^ was 
estimated following Wilson and Rood (|1994p . and assum- 
ing an [H2]/[C0] abundance ratio of 10^ (e.g. Scoville et 
al. 1986). The excitation temperature, Tg^, was estimated 
from the SMA-MRAM 30m ^^CO peak brightness temper- 
ature, Tb, of the averaged blueshifted and redshifted emis- 
sion of the outflows separately, assuming that TB=J{Tex)- 
J{Tbg), where the background temperature Tbg is 2.7 K. 
The Tex values are in the range 40-50 K. 

As seen in Table [H the values of Mout, Mout, P, E, 
and F of outflow A are about 4-5 times higher than those 
of outflow C. On the other hand, the size of the outflow 
C is larger and the dynamical timescales are comparable. 
Even if we add the redshifted SE feature to the calcula- 
tions of outflow C, the values estimated are still smaller 
than those of outflow A. This suggests that outflow C is 
less massive and powerful, probably because it is powered 
by a less massive and luminous source. In fact, following 
the correlations between the outflow parameters and the 
bolometric luminosity of the powering source, Lboi, found 
by Lopez-Sepulcre et al. (2009) for a sample of high-mass 
star-forming regions, the source powering outflow A should 
have a Lboi of a few lO'^-lO^ Lq, while the one driving out- 
flow C should have a L^oi of a few lO'^-lO'^L©. Furuya 
et al. (|2002p estimated a size of the lobes almost twice the 
size estimated with our observations. The fact that the syn- 
thesized beam of their observations was '~4'.'5, made them 
more sensitive to extended emission than ours. The Mout 
and Mout estimates obtained by Furuya et al. (j2002p for 
outflow A are three times smaller than ours, while for out- 
flow C their values are quite similar to ours. The dynamical 
timescales are comparable. 



4.3.3. Episodic outflow 

One of the striking features about the morphology of out- 
flow C is its clumpy structure. This is particularly visible 
in Fig. 1151 where we show the blueshifted and redshifted 
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emission maps averaged for high outflow velocities. As seen 
in this figure, the outflow lobes, in particular the redshifted 
one, show discrete knots. The knotty appearance of the 
outflow with multiple clumps aligned along the flow axis 
resembles that of the high-mass outflow HH 80-81 (Qiu & 
Zhang 2009 ), where a series of CO bullets have been dis- 
covered. These authors propose that the bullets have not 
been created in situ from entrained or swept-up ambient 
medium, but are the result of episodic, disk-mediated ac- 
cretion. For outflow C, the morphology could be explained 
in terms of an episodic outflow where the knots are made of 
swept-up ambient gas. In fact, in the position- velocity plot 
along the axis of this outflow, there is an increase in the 
extreme emission velocity at the position of all the knots 
for the lowest emission contour level. According to Arce &; 
Goodman (,2001b,l . the velocity structure at the position 
of the knots is characteristic of the prompt entrainment 
mechanism, which will produce the highest velocities at the 
position of the local CO maximum, and decreasing veloc- 
ity trend towards the source. The fact that the maximum 
outflow velocity does not increase linearly as a function 
of distance from the powering source (Fig. [T5)) indicates 
that there is a different "Hubble law" associated with each 
knot, corresponding with each mass-ejection episode. The 
extreme emission velocity for each mass-ejection episode 
will not necessarily be the same. This is because not all 
the outburst from an episodic source will necessarily have 
the same angle with respect to the plane of the sky or will 
accelerate the ambient gas to the same maximum velocity 
(Arce & Goodman F2001ap . 

5. Conclusions 

We analyzed millimeter data obtained with the SMA and 
the IRAM PdBI interferometers, and the IRAM 30-m tele- 
scope of the dust and gas emission towards the cores in 
the high-mass star-forming region G24. 78-1-0. 08. The syn- 
thesized beam of < 1" of the observations has allowed us 
to study with unprecedented high-angular resolution the 
structure and the velocity fleld of the cores, and the molec- 
ular outflows powered by the YSOs embedded in the HMCs. 

The millimeter continuum emission towards the cores 
Al and A2 has been resolved into additional cores. Core 
Al has been resolved into 3 cores named Al, Alb, and 
Ale, whereas core A2 has been resolved into 2 named A2 
and A2b. Interestingly, the cores are aligned in a southeast- 
northwest direction coincident with that of the molecular 
outflows detected in the region. This suggests a preferential 
direction for star formation in this region that could possi- 
bly depend on the direction of the magnetic field. It could 
also indicate star formation in filaments or sheets. The core 
associated with the HC Hii region is Al. The masses of the 
cores for which it has been possible to estimate the rota- 
tional temperature range from 7 to 22 Mq. The rotational 
temperatures range from 128 to 180 K. 

The velocity gradients towards cores Al and A2, first 
detected by Beltran et al. (PDUi|) in CHgCN (12-11), 
have been confirmed by the new SMA observations. 
Unfortunately, the sensitivity of the observations have not 
allowed us to better constrain the direction of the velocity 
gradient detected in CS (3-2) by Behran et al. (piM| . 

The CHgi^CN, i^CHaCN and C2H5CN observations 
have revealed the existence of 2 velocity components to- 
wards Al. The component at ^111 km s~^ is associated 



with the velocity gradient seen in CII3CN and peaks close 
to the position of the millimeter continuum peak and of 
the HC Hii region. The component at '^106 km s~^ peaks 
towards the SW of core Al and does not seem to be asso- 
ciated with a millimeter continuum emission peak. 

The two molecular outfiows in the region have been 
mapped in different tracers. The position-velocity plots 
along outflow A and the "'^^CO (2-1) averaged blueshifted 
and redshifted emission indicate that this outflow is driven 
by core A2. Regarding core Al, the position-velocity plots 
seem to indicate that this core is not powering any out- 
flow. However, due to the fact that the outflow structure 
is very complex towards cores Al and A2, we cannot to- 
tally discard the possibility that core Al could be pow- 
ering an additional outflow in the region. The outfiow C 
is highly coUimated and clumpy. The dumpiness and the 
^^CO (2-1) position- velocity plot suggest that this outflow 
could be episodic, and that the outflow clumps are made of 
swept-up ambient gas. 
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Fig. 4. Overlay of the emission averaged in the velocity 
interval [105, 115] km s'^of the CH3CN (12-11) X = 2, 
CH3CN (12-11) K = 6, CHa^CN (12-11) K=2, and 
^•^CHsCN (12-11) K = 2 transitions on the maps of con- 
tinuum obtained with the SMA and PdBI {bottom panel) 
towards G24. Contour levels are 3, 6, 9, 12, 15, 20, 30, and 
40 times a, where la is 15 mJy beam""'^ for all species ex- 
cept for "CH3CN (12-11) for which Icr is 8 mJy beam"!. 
The grayscale levels are the same as in Fig. [TJ The crosses 
mark the positions of the sources identified by Furuya et 
al. (|2002p and Beltran et al. (^004). The synthesized beam 
is shown in the lower left corner. 
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Fig. 9. Blueshifted {blue contours) and redshifted {red con- 
tours) i^CO (2-1) {top panel) and "CO (2-1) {bottom 
panel) averaged emission overlaid on the SMA VEX 1.3 mm 
continuum emission {grayscale) towards G24. The synthe- 
sized beam is 0'.'86 x 0'.'68 at P.A.= 54° for ^^CO and 
0'.'83 X 0'.'64 at P.A ~ 55° for ^^CO and is shown in the 
lower left corner. The emission has been averaged in the 
same velocity intervals as in Fig. [S] Contour levels are 3, 6, 
9, 12, 15, 20, and 30 times la, where la is 30 mJybeam"-'^ 
for ^^CO and ^^CO. Grayscale contour for the continuum 
emission is 3a, where la is 3.2 mJybeam^^. The crosses 
mark the positions of the sources identified by Furuya et 
al. (pnini) and Behran et al. 
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A2 in G24 with superimposed fit. Open circles and filled 
triangles correspond to the -'^^CHaCN (12-11) and the 
CH3CN (12-11) ws = 1 transitions, respectively. 
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Fig. 12. Overlay of the map of the 1.3 mm contmuum 
emission obtained with the very extended SMA configura- 
tion {contours) on that of the CH3CN (12-11) line velocity 
{colour scale). Contour levels are the same as in Fig. [1] 
Offsets are measured with respect to the phase center. The 
triangle marks the position of the HC Hii region in G24 
Al (Beltran et al. l2007p . The synthesized beam of the con- 
tinuum (line) emission is shown in the lower left (right) 
corner. 
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Fig. 13. {Top panels) ^^CHaCN (12-11) if =3 {black line), 
CHgi^CN (12-11) K=2 {dotted line), and C2H5CN (252,24- 
242,23) {gfO'V line) obtained by integrating the emission in- 
side the 3ct contour level area towards G24 Al {left) and A2 
{right). The vertical solid line indicates the T4.SR- {Bottom 
panels) Overlay of the averaged emission in the velocity in- 
terval (105, 107.5) km s'^ (component 1; left) and (108.75, 
114.25) km s-i (component 2; right) of the "CH3CN (12- 
11) if = 3 transition on the map of the 1.4 mm continuum 
emission obtained with the PdBI towards G24 Al and A2. 
Contour levels are the same as in Fig. [TJ The grayscale 
levels are 3, 6, 9, 12, 15, 20, and 30 times a, where la is 
3.7 mjy beam~^ for component 1, and 7.2 mJy beam~^ 
for component 2. The black triangle marks the position of 
the HC Hll region in G24 Al (Beltran et al. [2007|) . The 
synthesized beam is shown in the lower right corner. 
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Fig. 14. Overlay of the map of the 1.3 mm continuum, 
emission obtained with the very extended SMA configu- 
ration (contours) on the map of the first moment of the 
C2H5CN (252,24-242,23) second velocity component {colour 
scale). Contour levels are the same as in Fig.[TJ Offsets are 
measured with respect to the phase center. The triangle 
marks the position of the HC Hli region in G24 Al (Beltran 
et al. [2007p and the crosses the position of sources Al and 
Alb. 
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Fig. 15. {Top panel) Blueshifted {blue contours) and red- 
shifted {red contours) ^^CO (2-1) averaged emission to- 
wards G24 C. The emission has been averaged in the ve- 
locity intervals (58, 74.5) km s~^ {blueshifted emission) 
and (148, 167.5) km s~^ {redshifted emission). The molec- 
ular outflow map has been rotated —45°. The offsets are 
measured from the phase center, positive towards SE. 
Contour levels are 3, 4, 6, 9, and 12 times Icr, where Icr is 
15 mJybeam"-'^ (12 mJy) for the blueshifted (redshifted) 
emission. {Bottom panel) Same as bottom left panel of 
Fig. [TU] The vertical and horizontal solid lines indicate the 
position of core C and the Vi,sr, respectively. The position 
of each knot along the outflow axis is marked with a vertical 
dashed line. 
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Fig. 2. Upper and lower sideband SMA spectra obtained by integrating the emission inside the 3(t contour level area 
towards G24 A1+A2 (two upper panels) and G24 C (two lower panels) . The molecular lines are marked. The different K- 
components are marked with dashed lines in the upper (lower) part of each spectrum in the case of CH3CN (CHs^'^CN). 
The K, I quantum numbers of the vibrationally excited (fg = 1) CH3CN transitions arc indicated above the spectrum. 
Only the analyzed lines are labeled. 
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Fig. 5. Images of the emission averaged in the velocity interval (105, 115) km s^^for some observed molecular species 
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the sources identified by Furuya et al. (I2002p and Beltran et al. f|2004p . The synthesized beam is shown in the lower left 



Al 
,+ ,B, , , ,. 

_l I I I I L. 



16 



Beltran et al.: G24.78+0.08 




A 

a 

3 



Fig. 6. Combined SMA+IRAM 30-ni channel maps of the -'^^CO (2-1) emission towards G24. The contour levels range 
from 0.35 to 7.35 Jy beam~^ by steps of 1.4 Jy beam~^ and from 7.35 to 18.55 Jy beam~^ by steps of 2.8 Jy beam~^. 
The cross indicates the phase reference center of the SMA observations. The synthesized beam is shown in the lower 
right corner of the last panel. The central velocity of each velocity interval is indicated in the upper left corner of each 
panel. 
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Fig. 7. Same as Fig. [5] for the combined SMA+IRAM 30-m channel maps of the ^'^CO (2-1) hne. The contour levels 
range from 0.35 to 11.55 Jy beam"-'^ by steps of 1.4 Jy beam~^ 
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Fig. 8. Blueshifted (blue contours) and redshifted {red contours) ^^CO (2-1), "CO (2-1), C^^O (2-1), and SO (65-54) 
averaged emission overlaid on the SMA VEX 1.3 mm continuum emission (grayscale) towards G24. The maps have 
been reconstructed with a circular beam of 2" to enhance the extended emission. The emission has been averaged in 
the velocity intervals (82, 97) km s~^ {blueshifted emission) and (127, 142) km s~^ {redshifted emission) for ^^CO, (94, 
103) km s~^ {blueshifted emission) and (116, 127) km s~^ {redshifted emission) for ^'^CO, (99, 104.4) km s~^ {blueshifted 
emission) and (115.8, 122.4) km s^^ {redshifted emission) for C^^O, and (94, 104.2) km s~^ {blueshifted emission) and 
(115.6, 124.6) km s~^ {redshifted emission) for SO. Contour levels are 3, 6, 9, 12, 15, 20, 30, and 40 times la, where Icr 
is 80 mJybeam^^ for ^^CO, 35 mJybeam"^ for ^^CO, 18 mJybeam"^ for C^O, and 25 mJybeam"^ SO. Grayscale 
contour for the continuum emission is 3a, where la is 3.2 mJybeam"^. The spectral line synthesized beam is shown 
in the lower left corner. The crosses mark the positions of the sources identified by Furuya et al. (|2002p and Beltran et 
al. (.2004,) . 
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Fig. 10. Position- velocity plots along the direction with P.A.=:— 40° for outflow A {top and middle panels) and P.A.=— 45° 
for outflow C {bottom panels). The offsets are measured from the phase center, positive towards SE. Contour levels are 
3a, 6(7 9ct to 45cr by steps of 6ct where a is 0.033 Jybeam-i for ^^CO SMA, 0.33 Jy beam"! for ^^CO SMA+IRAM 30-m, 
0.027 Jybeam-i for "CO SMA, and 0.23 Jybeam"! for ^^CO SMA+IRAM 30-m. The SMA maps have the highest 
angular resolution. {Top and middle panels) The vertical solid and dashed lines indicate the position of cores Al and A2, 
respectively. The horizontal dashed line indicates the Vlsr- {Bottom panels) The vertical solid line indicates the position 
of core C. The horizontal dashed line indicates the Vlsr- 



